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The 15 N(p,7) le O reaction represents a break out reaction linking the first and second cycle of the CNO 
cycles redistributing the carbon and nitrogen abundances into the oxygen range. The reaction is dominated 
by two broad resonances at E p = 338 keV and 1028 keV and a Direct Capture contribution to the ground 
state of 16 0. Interference effects between these contributions in both the low energy region (E p < 338 keV) 
and in between the two resonances (338 < E p < 1028 keV) can dramatically effect the extrapolation to 
energies of astrophysical interest. To facilitate a reliable extrapolation the 15 N(p,7) le O reaction has been 
rcmeasured covering the energy range from E p =1800 keV down to 130 keV. The results have been analyzed 
in the framework of a multi-level R-matrix theory and a S(0) value of 39.6 keV b has been found. 
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I. INTRODUCTION 

The energy production and nucleosynthesis in stars is characterized by nuclear reaction sequences which determine 
the subsequent phases of stellar evolution. Energy production during the first hydrogen burning phase takes place 
jj^ ■ through the fusion of four protons into helium. This occurs either through the pp-chains or the CNO cycles. The pp- 
chains dominate hydrogen burning in first generation stars with primordial abundance distributions and in low mass, 
M<1.5 Mq, stars. The CNO cycles dominate the energy production in more massive, M>1.5 M Q , second or later 
generation stars with an appreciable abundance of CNO isotopes. The CNO cycles are characterized by sequences of 
radiative capture reactions and /3-decay processes as shown in Figure [l|. 

At stellar temperatures the 14 N(p,7) 15 reaction is the slowest process in the cycle, defining the time scale and the 
overall energy production rate [2|-|4j. This reaction is therefore of importance for the interpretation of CNO burning. 
The proton capture by 15 N is relevant since it is a branch point linking the first CNO or CN cycle with the second 
CNO, or NO cycle as shown in Figure [l[ This branching has always been a matter of debate since both reactions 
are characterized by strong low energy resonances. 

The reaction rate of 15 N(p,a) 12 C is determined by two broad low energy s-wave resonances at E p = 338 keV 
and 1028 keV, populating the J T = 1~ states at 12.44 and 13.09 MeV, respectively, in the 16 compound nucleus. 
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FIG. 1: Diagram of the CNO bi cycle [U 

There have been a number of low energy measurements [(| , Q , and [|| which provide the basis of the present rate 
in the literature [3, [13]. Recently, three lower energy points were derived from an indirect "Trojan Horse Method 
(THM)" [ll| which are consistent with low energy data j3-[3]. 

The competing 15 N(p,7) 16 reaction decays predominately to the ground state of 16 and exhibits the same two 
resonances but in addition is expected to have a strong non-resonant direct capture component [13]. The presently 
available low energy cross section data from [TU, [l4|, and [I3| differ substantially at lower energies. This poses 
difficulties for a reliable extrapolation of the cross section towards the stellar energy range. An extrapolation was 
performed fl2j ] using a two level Breit-Wigner formalism taking into account the direct capture contribution. The 
present reaction rate for 15 N(p,7) 16 [13] relies entirely on the predictions of [13]. 

To reduce the uncertainty in the strength of the direct capture term single-particle transfer reactions have been 
performed [l5[ to determine the proton Asymptotic Normalization Coefficient (ANC) for the ground state of 16 0. 
With this ANC value = (192 ±26) fm _1 , R-matrix fits to the 15 N(p,7) 16 data have been performed [H] which 

resulted in smaller values for the low energy cross section than those obtained by fl2j . This result was confirmed by 
an independent R-matrix analysis of the existing data [r|. This conclusion was furthermore supported by a study 
in which both the (p,a) and (p,7) reactions were analyzed simultaneously with a multi-level, multi-channel R-matrix 
formalism [l7j . 

Low energy data points were extracted from a re-analysis [lH of a study of 14 N(p,7) 15 performed at the LUNA 
underground accelerator facility at the Gran Sasso laboratory. This experiment was performed using a windowless 
differentially pumped gas target with natural nitrogen gas. For detecting the 7-ray signal, a large BGO scintillator 
detector was used to observe the characteristic 7 decay in summing mode [l9j |. Since the natural abundance of 15 N 
is low, the yield of the 15 N(p,7) signal was weak and overshadowed by beam induced background yield from proton 
capture reactions on target impurities and from the 14 N(p,7) reaction. The measurements were therefore limited to 
energies below 230 keV. The proposed cross section results are clearly lower than the values obtained by [13] but 
also slightly lower than the results of [13] and 14 1. Given the strong background conditions, the results could not be 
normalized to the known on-resonance yield at 338 keV and systematic errors in the data cannot be excluded. 

Because of the inconsistencies in the existing data and the uncertainties of an R-matrix analysis based on these 
existing data, we have performed a new study of the 15 N(p,7) 16 over a wide energy range using high resolution Ge 
detectors. In the following section the experimental approach at the Notre Dame Nuclear Science Laboratory (NSL) 
and the LUNA II facility at the Gran Sasso underground laboratory will be discussed. This will be followed by a 
discussion of the experimental data. In the last section the stellar nuclear reaction rate based on the present data will 
be calculated and compared with existing results. 
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II. EXPERIMENTAL SETUP 



A. Accelerators and Experimental Setup 



The experiment was performed at two separate facilities. At the University of Notre Dame the 4 MV KN Van de 
Graaff accelerator provided proton beams in an energy range of 700 to 1800 keV with beam intensities limited to 
<10 fj,A on target because of the high count rate in the Ge detector from the 15 N(p,ai7) 12 C reaction. The energy 
calibration of this machine was established to better than 1 keV using the well known 27 Al(p,7) 28 Si resonance at 
992 keV pi]. The 1 MV JN Van de Graaff accelerator at Notre Dame was used in the range of 285 keV to 700 keV 
with protons beams of 20 /xA. The energy of the this machine was calibrated using the well known 15 N(p,c*i7) 12 C 
resonance at 429 keV (20 1. 

The LUNA II facility 2JJ, located in the Gran Sasso National Laboratory, uses a high current 400 kV, Cockroft- 
Walton type accelerator. The accelerator provided proton beam currents on target of up to 200 fiA in the energy 
range of 130 to 400 keV. In addition to the high beam output, the accelerator is extremely stable, and the voltage is 
known with an accuracy of about 300 eV. 

The experimental setup in both experiments was very similar. The targets were water cooled and mounted at 45° 
with respect to the beam direction. At Notre Dame the position of the beam on the target was defined by a set of 
vertical and horizontal slits. The beam was swept horizontally and vertically across a target area of 1 cm 2 by steerers 
in order to dissipate power over a large target area. At LUNA the ion beam optics provided a de-focused beam on 
target and no beam sweeping was applied. To avoid the build-up of impurities on the target a Cu finger, cooled 
to LN2 temperatures, was placed along the inside of the beamline extending as close to the target as possible. In 
addition, a bias voltage of about -400 V was applied to the isolated cold finger to suppress the secondary electrons 
ejected from the target due to proton bombardment. 



B. Targets 



The Ti 15 N targets were fabricated at the Forschungszentrum Karlsruhe by reactive sputtering of Ti in a Nitrogen 
atmosphere enriched in 15 N to 99.95%. The stoichiometry was analyzed using Auger electron spectroscopy to confirm 
the composition. This test was performed on two target spots, one which had been exposed to beam and one which 
was not exposed. The results agreed within <2% with the nominal stoichiometry of 1:1. Isotopic abundances were 
experimentally verified by comparing the yield of the 14 N(p,7) 15 , 278 keV resonance [|[ from the enriched targets 
with that obtained using a target produced with a natural nitrogen gas. The results of this measurement showed an 
abundance of < 2% of 14 N for the thin targets corresponding to a 15 N enrichment of > 98% in agreement with the 
quote of the supplier. For the thick target used at LUNA, 14 N and 15 N enrichment of 17.4 ± 2.0% and 82.6 ± 2.0% 
were found, respectively, most likely caused by a contamination of the enriched gas during sputtering. 

The thicknesses of all TiN targets were measured using the narrow 15 N(p,ai7) 12 C resonance at 429 keV [2(J. The 
target used at Notre Dame had a thickness of 7.2 ± 0.3 keV at E p = 429 keV and the two targets used at LUNA had 
thicknesses of 9.5 ± 0.4 keV and 24.8 ± 0.5 keV at 429 keV, respectively. The stability of the Notre Dame targets 
was checked continuously during the course of the experiment using the 15 N(p,ai7) 12 C resonance. The thin LUNA 
target was monitored by rescanning the top of the 338 keV resonance in 15 N(p,7) 16 0. The thick LUNA target could 
also be monitored using the 14 N(p,7) 15 resonance at E p = 278 keV because of the large 14 N content (17%) of this 
target. Because of the relatively low power density delivered at Notre Dame, the target saw virtually no degradation 
over the experiment with an accumulated charge of 5 C (see Figure ^ and no yield corrections were necessary. During 
the LUNA experiment with significantly higher beam currents, the thickness of the thin target was reduced by 27% 
after an accumulated charge of 17 C and that of the thick target by 30% after an accumulated charge of 65 C (see 
Figure [3]). 



C. Detectors 



At the NSL, the 7-rays were observed using a HPGe Clover detector, which consists of four HPGe crystals contained 
in the same cryostat. This unique arrangement allows the separate detectors to be used in so called add-back mode [22[. 
At LUNA, a single crystal, 115% HPGe detector from Bochum, Germany was used for the detection of 7-rays. Several 
sample spectra are given in Figure |4j 

The primary advantage of the LUNA II facility is the low background environment created by the rock cover from 
the Gran Sasso mountains. The rock shields from cosmic rays and therefore decreases the 7- induced background from 
cosmic rays in the detector and has been shown [23[ to suppress the E 7 > 3.5 MeV background count rate in a Ge 
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FIG. 2: Notre Dame target (left) has a measured width of approximately 7.2 keV at E p = 430 keV. Squares, diamonds, and 
stars correspond to 0.006 C, 3.4 C, and 4.8 C of accumulated charge on the target, respectively. LUNA thin (right, diamonds) 
and thick (right, crosses) targets have measured widths of approximately 9.5 keV, and approximately 25 keV, respectively. The 
difference in stoichiometry between the two targets can be seen by comparing the plateau heights from the same resonance 
scan of each target (right). 



»? 0.20 



0.15 L 



g 

£ O.lOh 



^ 0.05 

5h o.oo 



I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 



"v -o. 





* * X 

♦°-x 


* 






*+ ■ 




/ * m 


* 








i i 1 i i i i 


A+ 

i i i i 1 1 Iii 


280 


290 300 310 



Proton Energy [keV] 



FIG. 3: Targets scans of thick target (18) using the 14 N(p,7) 15 0, E p = 278 keV resonance. The target was stable until about 
20 C, and went through significant deterioration. Squares represent the initial target scan, while down triangles represent the 
scan after 65 C on target 
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FIG. 4: Several sample spectra from the present 15 N(p,7) le O experiment. The spectra were taken with Notre Dame's KN 
accelerator (E p = 1300 keV, top), Notre Dame's JN accelerator (E p = 700 keV, middle), and the LUNA II accelerator (E p = 
130 keV, bottom). The 15 N(p,7) 16 ground state transition 7-ray peak is indicated in each spectrum with an arrow. 



detector by three orders of magnitude. In each experiment the detectors were set up at an angle of 45° with respect 
to the beam direction, allowing the position of the detector to be set as close as possible to the reaction position. The 
relative efficiency of the 7-ray detector systems was measured using radioactive sources along with well known capture 
7 reactions. At Notre Dame the relative efficiencies for high energy 7-rays were established using the 27 Al(p,7) 28 Si 
resonances at 992 keV QJ] and 1183 keV [H and the 23 Na(p,7) 24 Mg resonance at 1318 keV |||. The efficiency was 
extended to a 7-energy of 12.79 MeV using the n B(p,7) 12 C reaction at 675 keV and 1388 keV following the method of 
reference 25]. At LUNA the higher energy efficiencies were determined using the 278 keV resonance in 14 N(p,7) 15 
H] and the 163 keV resonance in 11 B(p,7) 12 C [2(j. The n B resonance has a small angular distribution [26[, and the 
correction for the relative intensity is less than 3%. 

While summing is not of concern for the ground state transition of the 15 N(p,7) 16 reaction itself, summing plays a 
significant role in the determination of the 7-efficiency. For this reason the efficiency measurements at both laboratories 
were carried out at several different detector-target distances and the data were simultaneously fitted for all distances 
following the procedure described by Imbriani et al. 0] (see Figure [5]). For the 1 cm and 5 cm measurements, there 
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y-ray energy [MeV] 

FIG. 5: Efficiency curves obtained at LUNA. Circles are measurements at a distance of 1 cm, squares at 3 cm, and diamonds 
at 20 cm. Closed symbols represent data without summing corrections, open symbols include summing corrections. 

is one point which has significant summing corrections. This point corresponds to the ground state transition in the 
E p = 278 keV 14 N(p,7) 15 reaction. While the ground state transition of this reaction has a small branching ratio, 
each of the other cascades have significant probabilities, which lead to strong summing effects [|[. 

III. EXPERIMENTAL RESULTS 

A. Cross Section Determination 

The excitation function for the ground state transition of the reaction 15 N(p,7) 16 has been measured in the energy 
range of 131 to 1800 keV. It consists of three distinct, overlapping sections, LUNA data from 131 keV to 400 keV, JN 
data from 285 keV to 700 keV, and KN data from 700 to 1800 keV. The experimental yield Y (number of reactions 
per projectile) at proton energy E p corresponds to the cross section er(E) integrated over the target thickness A: 

Y(E P )= ^ d E = a(E p ) T" (1) 

Je p -a £{E) Je p -a £{E) 

where e is the stopping power [27| , and f (E) represents the energy dependence of the cross section in the integration 
interval: 

a(E) = a(E p )f(E), (2) 

with cr(Ep) the cross section at E p . This reduces to the well known thin target yield equation if the cross section is 
constant over the target thickness, f(E)=l. However, at energies where the cross section varies significantly over the 
target thickness, the yield has to be corrected to extract the cross section. This correction factor is given by: 
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FIG. 6: Cross section for the ground state transition of the N(p,7o) O reaction versus center of mass energy. 
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and requires the knowledge of the energy dependence of the cross section in the energy interval A, f(E). This problem 
can be solved by an iterative method. In a first step f(E)=l is assumed resulting in an approximation of the cross 
section which in turn can be used to calculate a new correction factor. This process is continued until no change 
in the resulting cross section is obtained. This process quickly converges after 2 to 3 iterations. The resulting cross 
section is shown in Figure [5] 

An alternative method to determine the cross section was used in which the data were analyzed following the 
approach described in Q . Here, a brief review of the method is given, which consists of the analysis of the shape of 
the 7-line. 

The cross section, the stopping power of protons in TiN, the 7-ray efficiency, and the detector resolution all have 
an influence on the 7-ray line shape observed in a spectrum. The observed energy of the 7-ray is related to the proton 
energy by [J 
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The number of counts in channel i of the 7-spectrum, corresponding to the energy bin E 7 ^ to E 7 ^ + 5E 7 (<5E 
dispersion in units of keV per channel) is given by the expression 
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for E p i < E p (E p i = proton energy corresponding to channel i, E p = incident proton energy), where a{E P i) is the 
capture cross section, 77^(^71) is the 7-ray detection efficiency, e(E pi ) is the stopping power and is the branching 
of the associated decay. The conversion from E 7i to E pi includes the Doppler and recoil effects. The result is folded 
with the known detector resolution AE 7 to obtain the experimental line-shape. 

Finally, to infer the cross section in the energy window defined by the target thickness, the cross section was 
written as the sum of two resonance terms and a constant, non- resonant term as described in [28|. In these fits, the 
free parameters were the non-resonant astrophysical S factor and the gamma-ray background parameters. The results 
of both methods are in excellent agreement. 

The absolute cross section for the 15 N(p,7) 16 reaction has been measured on top of the two broad resonances at 
proton energies of 338 keV and 1028 keV using two independent well known reaction standards and two independent 
methods. All measurements were performed at three distances (d = 1 cm, 5 cm, 20 cm) to check for systematic errors. 
The exception to this protocol was the 15 N(p,7) 16 measurement at E p = 338 keV, where the yield at 20 cm was too 
low. In the first method the cross sections were determined relative to the thick target yield Y OQ ( 27 Al) of the well 
known 27 Al(p,7) 28 Si resonance at 992 keV: 



Foo( Al)= 2e(E R ) M ^ (6) 

with m (M) the mass of the projectile (target), A the DeBroglie resonance wave length and oj^r the resonance 
strength. Adopting the resonance strength of 1.93 ± 0.13 cV from [lj| [29|, H(| the cross sections values of er 338 =6.7 ± 
0.6 /ib and <7io28 =446 ± 40 /A are obtained. The errors include the statistical error (< 3 %), the error on the relative 
efficiency (2 %), the error on the relative charge measurement (2%), the uncertainty on stopping power values (5%), 
and the uncertainty of 7 % arising from the reference resonance strength of 27 A1. 

In the second method the cross sections were determined relative to the well known 429 keV resonance in 
15 N(p,ai7) 12 C (20[. In this case the cross section was calculated relative to the integral over the yield curve of 
the reference resonance, A re f . The resonance strength ujj is related to A re ^ by: 

"7=!^', (7) 

with n t the number of target atoms per cm 2 which can be calculated from the measured target thickness and the 
stopping power values. Using this relationship for A re ^ (Equation [7]), along with Equation [6] and the thin target yield 
approximation for the 15 N(p,7) 16 reaction, the following formula can be derived: 



X 2 (E R )\ (M 



X 



77(4.4MeVQ \ ( Y(E p ) \ 
V (12MeV) J \ A™f J W7i? ' 



with rj the relative 7-ray efficiency. In this approach the result is independent of all target uncertainties. The 
yield of the 4.4 MeV 7-rays from the 15 N(p,ai7) 12 C was corrected for the well known angular distribution [3TL l32j . 
Using the resonance strength of o>7,r=21.2 ± 1.4 eV from (33|, values of 0-338=6.4 ± 0.6 ^b and <Tio28=436 ± 44 fib 
are obtained. The errors include the statistical error (< 3 %), the error of the relative efficiency (2 %), the error of 
the relative charge measurement (2%), the error of A re f owing to the numerical integration and angular distribution 
correction (6%) and the uncertainty of 7 % of the reference value [33| . 

The weighted average of the results from the two methods gives final values of 0-338 = (6-5 ± 0.3) /ib (5%) and 
C1028 = (438 ± 16) ^b (4%) taking into account common error sources. The present results are compared to previous 
values I13-U4II in Table [T] and arc in very good agreement. An exception is the cross section at 338 keV from Rolfs and 
Rodney [12| which is about 50% higher than the other values. This discrepancy will be discussed in more detail in the 
next section. Data from LUNA and the Notre Dame JN accelerators were normalized to the low energy resonance, 
while data from the Notre Dame KN accelerator were normalized to the higher energy resonance. Good agreement 
between the data sets was found at the overlapping energies around 700 keV. 



B. Branching Ratios 

In previous experiments small cascade transitions have been observed for the broad resonances at 338 keV and 1028 
keV [2fJ]. For the 338 keV resonance a branching ratio of (1.2 ± 0.1) % for the transition to the 6050 keV level was 



9 





Present Rolfs 


Brochard 


Hebbard 




[12] 


M 


[13] 


unit:[keV] 


[fib] [/ib] 
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338 


6.5 ± 0.3 9.6 ± 1.3 


6.3* 


6.5 ± .3+ 


1028 


438 ± 16 420 ± 60 


490* 





TABLE I: Summary of present resonance cross sections in comparison to previous results. *No uncertainty is given. ' error 
derived from given reproducibility of 5%. 

found, in good agreement with the literature value of 1.2 ± 0.4 % [2(|. For the 1028 keV resonance, branches to the 
6050 keV level (1.0 ± 0.3) % and to the 7118 keV level (2.8 ± 0.3) % were measured, and were also in good agreement 
with literature value of (1.4 ± 0.4) % and (3.1 ± 0.8) %, respectively [Io|. 

C. Angular Distributions 

In previous analyses of the 15 N(p,7 ) 16 O reaction [12, EH > a direct capture component was required to fit the 
data. In each of these cases, only the 1.; = initial wave was included in the calculation. This component represents 
s-wave capture, and therefore interferes with the resonant component of the cross section. The s-wave components 
yield isotropic 7-ray distributions fl2l l28j. However, an li = 2 direct capture component is allowed for El transitions, 
and can contribute to the reaction, possibly introducing some anisotropy in the gamma decay. In looking for this 
signature, the set-up was tested using the 1 MeV resonance, which is reported to be isotropic [l2|. Measurements were 
made with the central Clover axis at several angles relative to the beam direction with a nominal distance of 5 cm. To 
extract additional data points, the add-back feature was used to treat the left and right halves of the Clover system 
as separate detectors, where Monte Carlo calculations using the code GEANT4 [34| yielded an effective angular offset 
of each half from the central detector axis of 15°. Measurements with a calibrated 137 Cs source, and the isotropic 278 
keV resonance of 14 N(p,7) 15 made corrections to the absorption effects possible. The results confirmed the isotropy 
of the angular distribution of the ground state decay 7-rays from the 1028 keV resonance in 15 N(p,7 ) 16 O. 

In the search for the 1/ = 2 component, angular distribution measurements of 15 N(p,7o) 16 were made at a proton 
energy of 540 keV at detection angles of 0°, 45°, and 90° relative to the beam axis. With two Clover segments, this 
resulted in 6 data points. Two of these points were excluded due to large absorption effects from the shape of the 
target chamber. The resulting angular distribution was fit with a function of the form W(8) cp P 2 (cos6>), 

where P2 is the L = 2 Legendre polynomial, and a^ p = Q2 &2, where Q2 is the usual geometrical correction factor 
due to the finite size of the detector [35|. The Q2 term can be directly determined using the 16 0(p,7) DC — > 495 keV 
reaction (28[. This reaction has a well known angular distribution of the form W{9) = sin 2 (6). Plotting the measured 
W(8) with respect to sin 2 (8), the deviation of the slope from 1.0 gives directly the Q2 value. Using the left and right 
Clover halves to measure the 16 0(p,7) 17 F DC^495 keV a value of Q 2 = 0.975±0.020 was found. This results in an 
a2 value for the 15 N(p,7 ) 16 O reaction at 540 keV of 0.08 ± 0.10, which is consistent with isotropy. 

IV. R-MATRIX ANALYSIS 

An R-matrix analysis was performed which mirrored the procedure of previous analyses [T2I [l5l [l6j . This analysis 
included the two broad 1~ resonance levels and a direct capture contribution. The analysis was performed with the 
multi-channel R-matrix code AZURE [36| . Details of the theory and the nomenclature are given in [36[ and references 
therein. The best fit results can be seen in Figure [7] where the cross sections have been converted to the astrophysical 
S-factor. This fit results in an S(0) value of 39.6 keV b. 

Looking closer at the low energy region (see Figure [5]) and comparing the results to previous experiments, the low 
energy data of [12] are inconsistent with the present results. While this data set agrees very well at higher energies it 
starts to deviate below a proton energy of 400 keV. There is no obvious reason for this inconsistency but it should be 
noted that the low energy data of Rolfs and Rodney [12] carry a significant uncertainty in this region. The data of 
Hebbard [l3[ are in reasonably good agreement above 230 keV. The data of Brochard [3] show good agreement except 
for the three data points below the 338 keV resonance which show a large scatter. The recent Bemmerer data [l8| are 
systematically lower than the present data, and also seem to have a slightly different energy dependence. 

The parameters for the best fit are given in Table [TTJ With the exception of the first resonance, most parameters 
are in fair agreement with previous results [TBI . [l6| . The width of the first resonance is dominated by the alpha 
width [20l | . This parameter is, therefore, well constrained by the present (p,7) data. The resonance strength, however, 



10 




0.5 1 1.5 

Center of Mass Energy [MeV] 

FIG. 7: Experimental S factor for the ground state transition of the 15 N(p,7o) le O reaction shown together with the best fit 
results of an R-matrix calculation using the code AZURE. 



Reference 


Ei 


Tpi 


7a 01 


771 ( int ) 


r 71 e 2 


2 2 
7p2 7ao2 


77 2 (int) 


1^72 


C^Ifm 1 / 2 ] 


Og.s 




r = r 


1 P = 


la = 1 




r = r 


l p = l a = 1 






s,d— >p 




Present 


12.438 


52.8 


13.5 


51.3 


33.8 13.087 


309.1 5.0 


34.1 


38.7 


23.22 


0.608 


[15] 


12.439 


280.9 


12.5 




8.8±1.5 13.089 


271.4 6.1 




50±8 


13.85 




[16, RR] 


12.452 


93.6 


13.5 


38.0 


13.111 


416.0 0.812 


0.5 






1.944 


[16, HH] 


12.447 


355.2 


10.6 


7.2 


13.087 


265.2 5.4 


56.2 






0.569 



TABLE II: Best fit R-matrix parameters for the 15 N(p,7) 16 reaction (E a in MeV, ^ in keV). Channel radii of a p = 5.030 fm 
and a a = 6.500 fm were used in the fit. The boundary condition BC in the fit was chosen to be equal to the shift function 
at the energy of the lowest J 7r = l~ resonance. Because of this choice of BC, the F 71 (in eV) can be calculated directly from 
the reduced width amplitude while T 72 is calculated from the Barker transformation of the reduced width amplitude [36| . The 
parameters shown in this table are mainly intended to enable the reader to reproduce the fit presented in this manuscript, and 
should not be regarded as final. With respect to the physical interpretation of these fit parameters, a detailed discussion will 
be presented in a forthcoming paper [37l j. 

is determined by the product of the proton and 7 width. Without including proton scattering data in the fit, these 
parameters are not well constrained. However, it should be noted that tests showed that this ambiguity has no 
influence on the extrapolation of the data. In these tests, the proton width (j Vl ) was fixed at different values while 
allowing the other parameters to vary. For the higher energy resonance the proton and alpha width are comparable, 
which provides more of a constraint on the parameters. In addition, the present result for the ground state ANC is 
significantly larger than the value of which cannot be attributed to the choice of radius (see below). While these 
differences do not have an impact on the extrapolation of the S-factor to lower energies, a more thorough analysis is 
warranted for the interpretation of the R-matrix parameters. This will be addressed in a forthcoming publication (37j 
where the results of simultaneous multi-channel fits to all relevant reaction channels will be presented. 
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FIG. 8: Low energy region of present data set (filled circles) shown with AZURE R-matrix calculation along with previous 
measurements of Rolfs [12] (open squares), Hebbard [l3|] (filled squares), Brochard [3|(open triangles), and Bemmerer [l8l] (open 
diamonds) in the region of - 400 keV 



In the present parameter space the only non-s-wave contribution arises from the d-wave component of the direct 
capture. Using the parameters for the best fit the 7-ray angular distribution was calculated at the minimum of the 
cross section between the two resonance at 540 keV yielding a value of a,2 = 0.034. The result is consistent with the 
experimental upper limit of 0.18 (see Section UTT1 C). 

The sensitivity of the best fit was tested against several key parameters. In multi-parameter fitting, uncertainties 
of specific parameters are determined in terms of confidence regions. For a nine parameter fit (E>,,7i for both levels 
plus ANC or 0g. s .) a 70% confidence region for one of the parameters is defined by the range where \ 2 < Xmin + 
10 [Hj]. The degrees of freedom is 102 (113 data points, 9 parameters) resulting in a reduced \ 2 °f 1-8. 



Radii of r„ = 5.03 fm, and r a = 6.5 fm, 



tested, and the results are given in Figure 9(a 



were taken from [161 ]. The dependence of the fit on the proton radius was 
Any proton radius value between 4 and 6 fm is considered acceptable. 



Variation of the radius over this range corresponds to only a 5% uncertainty in the S(0) extrapolation. This test also 
showed that the ANC is not very sensitive to the choice of the radius. The best fit gives a ground state ANC of (23 
± 3) fm -1 / 2 (corresponding to a reduced width amplitude of 8 g . s . = 0.61). The error associated with the ANC was 
determined by fixing the ANC at different values, and finding a new best fit. The results of this procedure can be 
seen in Figure 9(b) Even though the results give a 13% uncertainty for the ANC, the variation in S(0) from this 
procedure is only 4%. 

To evaluate the uncertainty of the fit itself, fits were performed which were forced to result in different S(0) values 
by inserting a fake data point at 1.5 keV with extremely small errors. The small error of this data point forces the 
calculation to match this fictitious data point and thus vary the extrapolation. By varying the fictitious point and 
observing the change in the x 2 ot , the experimental uncertainty of the extrapolation can be determined (see Figure ITU1) . 
This results in an uncertainty of ± 0.6 keV b. Including the 5% error of the absolute cross section give s a final S(0) 
value of (39.6 ± 2.6) keV b. This value is compared in Table HTT1 with pre vious extrapolations [H, M [l5|4l7|. The 
present result is in good agreement with the original value of Hebbard jl3l ] , the Barker [l6j analysis of the Hebbard 
data (labelled HH in Table IIII1 and the results of Mukhamedzhanov et al. []1| . The extrapolation of the Rolfs and 
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FIG. 9: Dependence of AZURE fits on key parameters used in the fits. Parameters were varied noting that for a 9 parameter 
fit, a 70% confidence region is defined for w \ 2 — Xmin + 10 Even though a lower \ 2 was found with a radius of about 
4.4 fm, 5.03 was chosen so as to be more comparable to Barker [l|| . ANC was found to be (23 ± 3) fm _1/ ' 2 . Fixing the ANC 
to the value found in [l^| gives a \ 2 thai is much higher than the current best fit value. 



Rodney data [l2[ (Rolfs and Rodney [l2[, Simpson [TtJ and Barker [T(|, labelled RR in Table HIT]) are all significantly 
higher because of the larger low energy cross sections used for their fit analysis. 



Analysis 


S(0) 7 (keV b) 


Hebbard 1960 [13] 


32 


Rolfs 1974 [12] 


64 ± 6 


Barker 2008 (RR) [16] 


« 50-55 


Barker 2008 (HH) [16] 


« 35 


Mukhamedzhanov [151 


36.0 ± 6.0 


Present 


39.6 ± 2.6 



TABLE III: Summary of Previous Results for extrapolated S(0) values 



V. REACTION RATE 

Using the results from the AZURE extrapolations, the reaction rates can be numerically determined using the 
formalism outlined in jioj 

N A {av) = 3.73- 10V"^-i (9) 

x / S(E)e(- 2 ^e(- lim5E WdE 
Jo 

with \i the reduced mass, T the temperature in GK, and rj the Sommerfeld parameter; the results are in cm 3 mole -1 
s . The above function was numerically integrated for temperatures between T = 0.01 - 10 GK. The present data 
covers an energy range of E = to 1.8 MeV, which validates the integration up to a temperature of T = 1 GK. For 
higher temperatures, the S-factor curve must be extended to higher energies beyond what the present experimental 
data cover. We followed the procedure of the NACRE compilation [13] which handles this difficulty by equating all 
higher energy S-factor values with the highest energy data point where the cross section varies only slowly with energy. 
The results are given in Table IIVI The present rate is at lower temperatures up to a factor two lower than previous 
rates reflecting the change in the low energy S-factor. 
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FIG. 10: Variation of x 2 with extrapolated S(0) values (for details, see text). 



Following the example of the NACRE compilations, the above results were fit using the following parametrization 

N A < av > = ai 10 9 T-hxp[a 2 T-^ - (T/a 3 ) 2 ] (10) 
[1 + a 4 T + a 5 T 2 ] +a 6 10 3 T^ 
exp{a 7 /T) + a$10 6 T-hxp(a 9 /T), 

where the best fit parameters can be found in Table Ivl 



VI. CONCLUSION 



The result of this work clearly demonstrates that there are significant uncertainties in the low energy cross section 
data of radiative capture processes of astrophysical relevance, despite many decades of low energy reaction studies. 
These uncertainties affect directly our understanding and interpretation of solar and stellar hydrogen burning phe- 
nomena. In this case the new results influence primarily the leakage rate from the CN to the ON cycle in stellar 
burning via the 15 N(p,7) 16 radiative capture process, which is reduced by a factor of two compared to the previous 
rate used traditionally in CNO nucleosynthesis simulations. In particular, the change in rate will modify the equilib- 
rium abundance of 16 0, which is correlated with the leakage rate of 15 N(p,7) 16 from the CN cycle and the rate of 
16 0(p,7) 17 F in the NO cycle. However a detailed study of the astrophysical impact of the present measurement goes 
beyond the aim of the present work, but should benefit from recent studies of low energy reaction rates [36| . 

The reliability of stellar reaction rates depends critically on the quality of the experimental cross section data. 
Direct measurements of the reaction cross sections at the Gamow range of stellar burning have been successful in 
only a few cases of reactions between light nuclei. Considering the anticipated count rates for the CNO radiative 
capture reactions we will continue to rely for most cases on the extrapolation of low energy measurements into the 
Gamow range. The present analysis clearly demonstrates that this requires a two-fold approach, pursuing the direct 
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T9 [GK] 


Nl((T«) Present 


Na(<?v) n acre 


N.4<<T'U)C/188 






[10] 


® 


0.010 


2.284E-21 


4.33E-21 


3.93E-21 


0.015 


1.377E-17 


2.66E-17 


2.34E-17 


0.020 


3.322E-15 


6.50E-15 


5.61E-15 


0.030 


3.215E-12 


6.41E-12 


5.38E-12 


0.040 


2.456E-10 


4.96E-10 


4.08E-10 


0.050 


5.378E-09 


1.09E-08 


8.90E-09 


0.060 


5.681E-08 


1.16E-07 


9.36E-08 


0.070 


3.747E-07 


7.64E-07 


6.14E-07 


0.080 


1.784E-06 


3.63E-06 


2.91E-06 


0.090 


6.699E-06 


1.36E-05 


1.08E-05 


0.100 


2.104E-05 


4.23E-05 


3.37E-05 


0.150 


1.280E-03 


2.46E-03 


1.92E-03 


0.200 


1.908E-02 


3.52E-02 


3.09E-02 


0.300 


6.143E-01 


1.05E+00 


1.57E+00 


0.400 


4.682E+00 


7.48E+00 


1.45E+01 


0.500 


1.668E+01 


2.58E+01 


5.42E+01 


0.600 


3.908E+01 


5.90E+01 


1.27E+02 


0.700 


7.213E+01 


1.07E+02 


2.25E+02 


0.800 


1.165E+02 


1.68E+02 


3.39E+02 


0.900 


1.763E+02 


2.46E+02 


4.67E+02 


1.000 


2.599E+02 


3.50E+02 


6.18E+02 


1.500 


1.485E+03 


1.67E+03 


2.33E+03 


2.000 


4.735E+03 


5.05E+03 


6.61E+03 


3.000 


1.489E+04 


1.55E+04 


1.97E+04 


4.000 


2.428E+04 


2.27E+04 


3.12E+04 


5.000 


3.064E+04 


2.61E+04 


3.83E+04 


6.000 


3.435E+04 


2.86E+04 


4.19E+04 


7.000 


3.624E+04 


3.06E+04 


4.31E+04 


8.000 


3.699E+04 


3.23E+04 


4.29E+04 


9.000 


3.706E+04 


3.37E+04 


4.19E+04 


10.000 


3.674E+04 


3.51E+04 


4.04E+04 



TABLE IV: Table of reaction rates for 15 N(p,7 ) 16 O. Present rates, along with the NACRE [T3], and CA88 @] results are given. 



ai = 


0.523 


a4 = 


6.339 a 7 


= -2.913 


a2 = - 


15.240 


a 5 = 


-2.164 a 8 


= 3.048 


a3 = 


0.866 


ae = 


0.738 a 9 


= -9.884 



TABLE V: Best fit parameters for the 15 N(p,7o) le O reaction rate using the NACRE fitting formulations. 

reaction measurements to the lowest possible energies in a background shielded environment but also expanding 
the experimental range of the measurements to determine unambiguously the various reaction components of the 
radiative capture process. The latter step is essential for minimizing the uncertainties in the R-matrix analysis of the 
cross section and can be complemented by independent studies which explore independently the strength of specific 
"hidden" reaction components such as the direct capture through ANC measurements and analysis. 

The approach taken here for the study of the 15 N(p,7) 16 reaction has succeeded in combining both the efforts of 
improving on the extent and quality of the low energy cross section data in underground accelerator experiments. At 
the same time the study has improved on the detailed measurement of higher energy data providing a better constraint 
on determining the external capture component and its impact on the low energy extrapolation of the reaction cross 
section. The combination of these two complementary measurements successfully reduced the overall uncertainty in 
the 15 N(p,7) 16 reaction rate. 
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